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Abstract: In the framework of Soft-Collinear Effective Theory, the jet quenching pa- 
rameter, q, has been evaluated by adding the effect of Glauber gluon interactions to the 
propagation of a highly-energetic collinear parton in a medium. The result, which holds 
in covariant gauges, has been expressed in terms of the expectation value of two Wilson 
lines stretching along the direction of the four-momentum of the parton. In this paper, 
we show how that expression can be generalized to an arbitrary gauge by the addition of 
transverse Wilson lines. The transverse Wilson lines are explicitly computed by resumming 
interactions of the parton with Glauber gluons that appear only in non-covariant gauges. 
As an application of our result, we discuss the contribution to q coming from transverse 
momenta of order g 2 T in a medium that is a weakly-coupled quark-gluon plasma. 
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1 Introduction 

There is accumulating evidence that a new state of matter, called the quark-gluon plasma, 
is formed in heavy-ion collision experiments at RHIC and LHC. The properties of this new 
state can be studied by using different observables and phenomena. One of them is the 
so-called jet quenching. 

Jet quenching is the process by which a highly-energetic jet loses energy while traver- 
sing a medium. Jets are created in the early stages of the collision before the plasma has 
been formed. When the parton that fragments into the jet interacts with the medium, it 
may lose energy by different processes, such as gluon bremsstrahlung and pair production 
[1-10] (in [11] there is a comparison of different formalisms and further references). By 
comparing suitable jet observables from heavy- ion collisions with observables from p-p 
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collisions, it is then possible to gain information about the properties of the medium. Jet 
quenching has been measured by the CMS [12] and ATLAS collaborations [13] at LHC, and 
before by the PHENIX [14] and STAR collaborations [15] at RHIC. The measurements at 
RHIC have been analyzed in different theoretical models (see [16] for a review). A common 
feature of the models is that they need to consider the effect of the so-called transverse 
momentum broadening. 

Transverse momentum broadening refers to a process that happens when the parton 
forming the jet interacts with constituents of the medium that have a much lower energy 
than the parton. The interaction does not change (in first approximation) the energy 
of the parton but it changes the momentum component of the parton that is transverse 
to the initial jet direction. This process is commonly characterized by the so-called jet 
quenching parameter, q [3]. Theoretically, q can be related in a model- independent way 
to the expectation value of two Wilson lines oriented along one of the light-cone directions 
[2, 4, 17, 18]. A way to derive this result is by factorizing the high-energy physics of the 
jet from the low-energy physics of the medium. A convenient set up is provided by the 
Soft-Collinear Effective Theory (SCET), which is an effective field theory (EFT) suited to 
describe processes involving nearly massless highly-energetic particles like jets [19-24]. 

The expression of q in terms of Wilson lines oriented along the light-cone holds only in 
covariant gauges and is not, in general, gauge invariant. This can best be seen by choosing 
the light-cone gauge A + (x) = 0, which sets all Wilson lines equal to one. In [25], a similar 
problem was studied for semi-inclusive deep inelastic scattering (SIDIS); the result was 
extended by analogy to the jet quenching parameter case. The expressions obtained in the 
SIDIS case are manifestly gauge invariant, while the expression for q is gauge invariant 
only for some choice of regularization of the light-cone gauge singularity. In other studies, 
the expressions found for q are not gauge invariant; a gauge invariant expression for q can 
be found in [26], but limited to second order in the opacity expansion. 

In the last few years, some studies have addressed SCET in the light-cone gauge. These 
studies have stressed the relevance of transverse Wilson lines to restore gauge invariance 
in some observables and in the SCET Lagrangian [27-29]. We will see that, as conjectured 
in [18, 30], transverse Wilson lines will also be crucial to make the expression of q gauge 
invariant. Part of the difficulties in obtaining a fully gauge- invariant expression for q is 
related to the fields in the Wilson lines being path ordered but not time ordered [18]. Hence, 
fields at equal times located in different points of the Wilson lines are not contiguous, a 
fact that requires some care when dealing with gauge transformations. 

The main aim of the paper is to derive a gauge invariant expression for q. This will 
be achieved by explicitly calculating in SCET and in light-cone gauge the scattering of a 
collinear highly-energetic parton on a background of Glauber gluons. Glauber gluons are 
gluons that do not modify the collinear nature of the parton, however, they can signifi- 
cantly change the transverse component of its momentum. Glauber gluons have first been 
included in the SCET Lagrangian in [31]; as it will turn out, the form of the leading-power 
Glauber interaction will depend on the gauge used. Among others, Glauber gluons are 
responsible for the transverse momentum broadening of the jet. Other responsibles are 
soft and collinear gluons [18]. We will argue that, at lowest order, by resumming Glauber 
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gluons the possible effect of soft gluons is automatically taken into account. Collinear 
gluons provide an independent source of transverse momentum broadening (and energy 
loss). In this work, we will not consider collinear gluons and, more in general, the effect 
of fragmentation into collinear partons. The fact that resumming Glauber and soft gluons 
gives rise to a well-defined gauge-invariant expression appears to justify, at least to the 
order at which we are working, separating their effect from the one of collinear gluons. 

Finally, having obtained a gauge invariant expression for q will allow us to make use 
of the arguments developed in [32, 33] and relate, under some assumptions, q to thermal 
field averages evaluated numerically in lattice gauge theories. In particular, in the case of 
a weakly-coupled quark-gluon plasma of temperature T, we will discuss the contribution 
to q coming from transverse momenta of order g 2 T. 

The paper is organized in the following way. In section 2, we set the general frame- 
work of the computation. Specifically, we will discuss the power counting of the SCET 
Lagrangian in different gauges, as well as the properties of the gauge fields at light-cone 
infinity. In section 3, we define and compute the transverse momentum broadening of a jet 
and the jet quenching parameter. We shortly review the calculation in covariant gauges, 
while we detail the new calculation in the light-cone gauge A + (x) = 0. Finally, we present 
an expression valid for arbitrary gauges. As an application of our result to the case of a 
weakly-coupled plasma, in section 4, we use the gauge invariant definition of q to discuss 
the contribution from transverse momenta of order g 2 T. In section 5, we conclude. 

2 Effective field theory, Feynman rules and power counting 

In this paper, we aim at obtaining a gauge invariant expression for the jet quenching 
parameter. Our work has been inspired by [18], which, in turn, follows from the study 
of [31]. In these references, SCET was extended to include so-called Glauber gluons, 
which will be defined in the next section, and used as the framework for the calculation 
of jet broadening. This is advantageous, because in an EFT framework quantities may 
be calculated by systematically expanding in the different relevant kinematical regimes 
already at the Lagrangian level and with a definite power counting. Effective field theories 
also offer the advantage that some classes of potentially large perturbative contributions 
may be resummed via renormalization group equations. Here, we will make a limited use 
of the potentialities of EFTs as we will focus on few kinematical regimes and not solve 
renormalization group equations. Still we believe that SCET provides a transparent power 
counting and a suitable framework for further improvements. 

Differently from the above references, we will assume a general gauge scenario and 
explicitly work out the light-cone gauge case. This will require a different power counting, 
the introduction of new Feynman rules and a particular care in treating gauge fields at 
infinity in one light-cone direction. We will discuss these issues in the rest of the section. 

2.1 Energy scales and degrees of freedom 

We are looking at highly-energetic jets propagating in a medium made of low-energy parti- 
cles. Hence, such systems are characterized by a large and at least one small energy scale. 
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The large scale is the momentum Q of the primary jet particle (quark or gluon) in one 
light-cone direction (or the energy, which is of the same order). The smaller scale is the 
temperature, T, or any energy scale characterizing the medium in non-thermal systems. 
This gives rise to a small dimensionless ratio A, e.g. A = T/Q <C 1. 

We will consider a parton moving along the light-cone direction n with initial momen- 
tum qo = (OjQjO). 1 Furthermore, we require the parton to travel through the medium 
without going far off-shell or losing much energy through radiation; more specifically, the 
virtuality of the parton in the final state should be of the order of Q 2 X 2 or smaller. Since 
we are looking at transverse momentum broadening through interaction with the medium, 
we further require the parton to acquire a large transverse momentum of the order of QX. 
In short, we will consider the propagation through the medium of a parton whose final mo- 
mentum scales like Q(X 2 , 1, A); such a parton is called collinear. The parton will eventually 
fragment into a jet, corresponding to a narrow cone of particles with a large energy and a 
much smaller invariant mass. 

We may classify particles interacting with the parton according to the size of their 
momenta. We call hard all particles whose momentum square, i.e. virtuality, is of or- 
der Q 2 , soft, particles whose momenta scale like Q(X,X, A), Glauber, particles whose mo- 
menta scale like Q(A 2 ,A 2 ,A) or Q(X 2 ,X,X) and ultrasoft, particles whose momenta scale 
like Q(X 2 , A 2 , A 2 ). 2 Therefore, we assume that particles in the medium have a typical soft 
momentum. 

Energy loss through radiation can be induced by the emission of collinear partons. This 
is the primary ingredient for the calculation of jet fragmentation. Following analogous pre- 
vious analyses, we will not consider here the effect due to the coupling of collinear particles, 
keeping however in mind that such an effect is not suppressed by any power counting. A 
first study including collinear radiation effects on the propagation of an energetic parton 
in a medium can be found in [35]. 

If coupled to a collinear particle, soft modes give rise to a so-called hard- collinear 
particle with momentum Q(X, 1, A), which is off shell by a momentum square of order Q 2 X. 
It has been argued in [18] that because the virtuality of this particle is larger than the 
one of a collinear particle, which is Q 2 X 2 , the effect due to the interaction with soft modes 
is suppressed in the strong coupling constant. In this work, we will, however, explicitly 
include soft and hard-collinear modes and show that their possible effect can be cast in the 
same gauge-invariant expression that encompasses the effect of Glauber gluons. 

The most relevant contribution to the transverse momentum broadening of a single 
collinear parton originates from the interaction with Glauber gluons. In the following, we 
will consider their contribution together with the one of soft gluons. 

Our notation is q — (q + , q~ ,q±), with q + = n ■ q, q~ — n ■ q, and the light-cone directions are given by 
fi = ^75(1, 0, 0, —1), n — -^75(1, 0, 0, 1), such that n • n = 1. The transverse momentum component is q±. 

2 The original definition of Glauber particles includes only particles with momentum scaling like 
Q(A 2 ,A 2 ,A). Following [26], we extend our definition to include also particles with momentum scaling 
like Q(X 2 , A, A) whose importance has been recently stressed in [34]. 
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2.2 SCET 

We restrict our analysis to the case of a highly-energetic primary parton that is a light 
quark. The SCET Lagrangian that describes the propagation of a massless quark in the 
light-cone direction n is 

Cn = Iniin ■ D i n + £ n ity ± . 1 ityx. iin , (2.1) 

2in ■ V 

where denotes the quark field and iD^ = id^ + gA^. The Lagrangian (2.1) is suited to 
describe the propagation of a collinear quark interacting with soft, collinear, Glauber or 
ultrasoft gluons. Hence, the quark field £ n is collinear when acting on the initial and final 
states, but may describe both hard-collinear and collinear virtual quarks. 

We now proceed through the following two steps. First, we rescale the collinear field 
£n by the large momentum component Q, according to ^ n — > e~ l ® x This implies that 
the residual momentum of the quark along the light-cone direction n, which is due to the 
interaction with the medium, is now of order QX or smaller. Second, we exclude from the 
Lagrangian collinear gluons. As we argued in the previous section, this last requirement 
does not rely on the power counting, but on the convenience to split the calculation in a part 
that deals with one collinear particle in the final state and a part that deals with more than 
one collinear particle in the final state, while postponing the latter for future considerations. 
Under the above conditions, we have 2in ■ D = 2Q + residual momenta of order QX or 
smaller. This implies that the last operator in (2.1) may be expanded into local operators, 
which leads to 3 

D 2 - gG^ u 

Cn = Zniin- D €n + €n-^ Hn + in i ^Q~ 1^ Hn + • • ■ , (2.2) 

where the dots stand for higher-order terms in the A expansion. The SCET Lagrangian 
(2.2) amounts to having integrated out from QCD hard modes at leading order in the 
strong coupling constant, as well as the small components of the collinear fields. 

The SCET Lagrangian that describes the propagation of a collinear particle and its 
interaction with Glauber gluons has been introduced in the context of jet scattering on cold 
nuclear matter in [31]; Glauber gluons have also been found necessary for the consistency 
of factorization proofs in certain processes like Drell-Yan [36]. 

2.3 Power counting and leading order Lagrangian in covariant gauge 

In the SCET Lagrangian, collinear and hard-collinear quark fields, £ n ( x )i scale in the same 
way. This follows from the fact that both the inverse propagator in momentum space of a 
hard-collinear quark and the typical four-momentum region occupied by a hard-collinear 
quark are enhanced by the same factor, A, with respect to the collinear quark case; hence, 
hard-collinear and collinear quark propagators scale in the same manner in position space 



We have used 

l/) ± ]p ± = -D ± ■ D ± - \gG^ liilv , 



where G x " = i[D^_, D±]/g is the gluon field strength. Here and in the following, v± denotes a vector such 
that u x = (0,v\v 2 ,0) and v\ = (v 1 ) 2 + (v 2 ) 2 = -<«x M - 
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and so do hard-collinear and collinear quark fields. The operators n ■ d and = d±i 
scale like QX 2 and QX respectively when acting on a collinear field Cn{x), and both scale 
like QX when acting on a hard-collinear field £n(x). Soft-gluon fields scale like QX and 
ultrasoft-gluon fields scale like QX 2 , for they are homogeneous in the soft and ultrasoft 
scale respectively. In contrast, the power counting of Glauber gluons depends on the 
gauge. The equations of motion require A + (x) to scale like n ■ d when acting on a collinear 
quark, or smaller. We assume A + (x) ~ QX 2 . In a covariant gauge, if the gluon field is 
coupled to a homogeneous soft source, this also implies A±(x) ~ QX 2 . In [26, 31] other 
sources have been considered that lead to different power countings. 4 

According to the power counting, the leading-order SCET Lagrangian in a covariant 
gauge is 

V 2 

£-n '° OV = iniin- D^n+^n-^i^n- (2.3) 

The first term on the right-hand side may either involve the interaction of two collinear 
quarks and a Glauber gluon or the interaction of a collinear quark with a soft gluon and 

a hard-collinear quark. Both interactions are of order one: J d 4 x £s A^ (x) £^ ~ 1. This 

implies that a collinear quark may exchange in the medium an arbitrary number of Glauber 
or soft gluons, all contributing to the same order in A to the scattering amplitude. The 
second term on the right-hand side of eq. (2.3), instead, involves only collinear quarks, 
the hard-collinear terms being suppressed by one power of A. From the leading-order 
Lagrangian, it follows that only the propagator of a collinear quark, of a hard-collinear 
quark and the interaction vertex with the A + component of a gluon matter. The propagator 
of a collinear quark carrying momentum q is 

the propagator of a hard-collinear quark carrying momentum q is 

(2.5) 



2q+ + it 



4 Soft sources have been considered both in [18] and [26]. In [18], for Glauber gluon momenta scaling 
like Q(A 2 , A 2 , A) they find that all components of the Glauber field scale like QA 2 , while in [26] for Glauber 
gluon momenta scaling like Q(A 2 ,A,A) they find that all components of the Glauber field scale like QA. 

Both results, however, seem to follow from an incorrect handling of the phase-space integral J d 4 pS(p° — E) 

where E generally scales like QA, which is the typical energy of the soft source. The integral, correctly 
evaluated, vanishes over the momentum region p ~ Q(A 2 ,A 2 ,A), for p° ~ QA 2 < £ ~ QA, while it is 

j d 4 pS(p° — E) — J dp + dp~ d 2 p± 8(p + /V2 + p~ /V2 — E) — \/2 J dp + d 2 p± ~ Q 3 A 4 over the momentum 

region p ~ Q(X 2 , X, A). Hence it is only Glauber gluons whose momenta scale like Q(A 2 , A, A) that couple at 
leading power to a soft source. Using the above integral, it follows that all components of Glauber fields of 

this kind scale like QA 2 , which is the scaling adopted in this paper. The identification of J d 4 pS(p° — E) 

with J d 3 p is correct only for Glauber gluons whose momentum scales like Q(A 2 , A, A) and after replacing 
p 3 by the small momentum y/2p + . A similar observation can be found in [31]. 
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whereas the vertex of the n fJj A +a component of a gluon with a collinear or hard-collinear 
quark is 

igT a f^i. (2.6) 

The Lagrangian (2.3) contains, in principle, also ultrasoft gluons. It has been shown, 
however, that ultrasoft gluons decouple at lowest order from collinear quarks trough the 



field redefinition £n(x) — > P exp 



ig dyn- A us (x + , y, x±) 



£n(a;), where P stands for 



the path ordering operator and A us for an ultrasoft gluon field [22]. The field redefinition 
works for ultrasoft gluons because their transverse momentum is suppressed with respect 
to the transverse momentum of collinear quarks, so that, at lowest order, the kinetic 
energy operator, Vj_/(2Q), commutes with ultrasoft gluons. Note that, for the opposite 
reason (the transverse momentum of Glauber gluons is of the same order as the transverse 
momentum of collinear quarks), the field redefinition would not decouple, even at lowest 
order, Glauber gluons from collinear quarks. 

We conclude this section by observing that in a covariant gauge the free gluon propa- 
gator of a Glauber gluon may be approximated by 



d 4 xe lk - x {0\T(A^ l (x)A u (0))\0) ^ 

k i 



g^ u - ex — 



(2.7) 



where a is a gauge parameter. Equation (2.7) follows from k 2 ~ k\ ~ Q 2 X 2 and k + ~ 
QX 2 <C k± ~ QX. Hence, in the case of Glauber gluons and at lowest order in A, the gluon 
propagator does not depend on the light-cone momentum component k + . This also holds 
for background Glauber fields coupled to soft sources. 

2.4 Light-cone gauge 

The Lagrangian (2.3) cannot be the lowest-order Lagrangian in the light-cone gauge, 
A + (x) = 0, because, otherwise, the collinear quarks would decouple from the gluons and 
propagate freely. It turns out, indeed, that in light-cone gauge, Glauber gluons scale dif- 
ferently than in covariant gauges and, because of this, new interaction vertices between 
collinear quarks and gluons must be considered at lowest order. 

The different scaling of the gluon fields in the light-cone gauge A + (x) = is due to the 
fact that, in that gauge, the gluon propagator is singular for k + — > 0. As a consequence, 
A±(x) does not vanish for x~ = ±oo, like in covariant gauges, and may be conveniently 
decomposed into [27, 28] 

A\{x) = A c ° v '\x) + Ai n '\x) , (2.8) 

with 

A s i n > l (x) = 8(x-)A\(x + ,oo,x ± ) + 9(-x-)A\(x + , -oo,x ± ) . (2.9) 

The field A c ^ v ' i (x) contributes to the non-singular part of the propagator and vanishes at 
x~ = ±oo, while the field A^^lx) contributes to the singular part of the propagator. Note 
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that the theta functions, 9(±x~), are singular in momentum space for k + — > 0. They 
multiply (non- vanishing) fields evaluated at x~ = ±00. 

Before turning to the SCET Lagrangian, it is useful to consider some general properties 
of gauge fields at x~ = ±00 [37] . The energy of the gauge field is proportional to the integral 
over space of E 2 + B 2 , where E and B are the chromoelectric and chromomagnetic fields 
respectively. The energy is finite if both fields vanish at infinity. This, in turn, requires 
the fields to be pure gauge at infinity; so that at infinity we may write 

A fl (x) = -d^(x)n\x), (2.10) 
9 

where £l(x) is an SU(3) gauge transformation. For an infinitesimal gauge transformation 
il(x) ~ 1 — igcp(x), this implies 

A^x) = 8^{x) . (2.11) 

In the light-cone gauge A + (x) = 0, the condition at infinity (2.11) translates into the 
following condition for the transverse components of the gauge fields at x~ = ±00: 

A i ± (x + ,±oo,x ± ) = -V i ± (p ± (x + ,x ± ) . (2.12) 

The fields =t , solution of eq. (2.12), are 

(f> ± (x + ,x±) = - / dsl±- A ± (x + , ±00, x ± + l ± s) , (2.13) 
J— 00 

where, for convenience, we have chosen to integrate over a straight line going from x± —ool± 
to x±, l± being an arbitrary vector in the transverse plane. The solution does not depend 
on the chosen integration path: for a path variation can be written as a surface integral 
over a field strength tensor [38], which vanishes at x~ = ±00. The solution does depend 
instead on the starting point of the integration path, which, in our case, has been chosen 
to be x± — ool±. In a covariant gauge, because there is no preferred direction in space, 
this cannot happen, which is another way to see why the field has to vanish at x~ = ±00. 
In light-cone gauge, different starting points lead to different fields 4> that are related 
by a gauge transformation. Gauge invariance then requires physical observables, like the 
jet quenching parameter, to be independent of the choice of the starting point. We will 
explicitly show that this is indeed the case. 

2.5 Power counting and leading order Lagrangian in light-cone gauge 

In the light-cone gauge A + (x) = and for soft sources, the Glauber fields A~(x) and 
Af v (x) scale as in the covariant gauge case, i.e. like QX 2 . However, the singular part of 
the Glauber field A±(x) is enhanced by one power of 1/A, so that we have A s ^ n,1 (x) ~ QX. 5 
Soft and ultrasoft gluons scale like in the covariant gauge case and do not contribute at 

5 This may be traced back to the factor k I L /k + appearing in the Fourier transform of the ghion propagator 
(0| T (AT(x)A Bia -(0)) |0>; k\/k + is of order 1/A for Glauber gluons. 



-8- 



lowest-order to the Lagrangian Cj%. From eq. (2.2) and according to the above power 
counting, the leading-order SCET Lagrangian in light-cone gauge is then 

C-n =£nl$Th-o£n + £n 2 q Hn ■ (2.14) 

The gluons in the Lagrangian are Glauber gluons and the quark fields are collinear quark 
fields. Note that the operator proportional to does not appear at leading order because 
the fields (jx^ do not contribute to the field strength tensor, which, in turn, reflects the fact 
that the field strength tensor vanishes at infinity. 

According to the leading-order Lagrangian (2.14), in light-cone gauge, collinear quarks 
couple to (transverse) gluons in two possible ways: via a 1-gluon and a 2-gluon vertex. 
Keeping Glauber gluons as background fields, the Feynman rule for the 1-gluon vertex 
reads 

* , q' ± • gA™ (q' - q) + g^f^ - q) • q± 

--% — ft, (2.15) 



whereas the Feynman rule for the 2-gluon vertex is 

i i" if dV 




2QJ (2ttY 



gA s l ni (q" - q) gA s l nt (q -q)$. (2.16) 



3 Jet broadening 

The jet quenching parameter, q, can be defined as the average square transverse momentum 
with respect to the original direction of motion that a highly-energetic parton picks up 
while travelling a large distance L through the medium. Following the notation and the 
derivation of [18], q can be written as 

1 n2 \ _ 1 f d 2 k ± 2 



Q = = 1 J j^p k ± p ( k ±) , (3-1) 

where P(k±) is the probability for the hard parton to pick up a certain transverse momen- 
tum k±. The calculation of P(k±) and q proceeds then through the following steps. 

(1) First, we consider a highly-energetic parton of momentum qo = (0, Q, 0) in an 
initial state, |in), propagating in a box of length L (volume L 3 ) and interacting with an 
arbitrary number of background gluons from the medium leading to a final state, \k, a), 
made of a collinear parton of momentum k and polarization a. Squaring the amplitude 
and integrating/summing over the final state momenta/polarizations and colors, we obtain 
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where the initial state has been normalized relativistically (the initial state energy is Q/y/2) 
in a box of size L. T is the interaction part of the S-matrix. It is useful to write T = Y2 n T n , 
n being the number of gluon lines attached to the parton, and define 

d ■*A.'jfi 77 1 \ > f dk 



dk 2 



\/2QL 



3 £ 



(2vr)2|k| 



(in| \k,a)(k,a\ T n |in) 



so that 



d 2 k\ d 2 A r , 



dki 



(3.3) 



(3.4) 



(2n) 2 

A convenient way to calculate A is via the optical theorem. This amounts to computing 
twice the imaginary part of all diagrams made of an incoming and outcoming parton of 
momentum go interacting with an arbitrary number of background gluons. The quantity 
A mn is then the cutting diagram made of n gluons before the cut and m gluons after the 
cut, see fig. 1. In the quark case, cutting the collinear quark propagator amounts to picking 
up the discontinuity of eq. (2.4), while hard-collinear propagators do not contribute to the 
cut because they do not correspond to on-shell modes. Vertices and propagators on the 
left-hand side of the cut are those defined in sections 2.3 or 2.5; vertices and propagators 
on the right-hand side of the cut are the complex conjugated of those on the left-hand side. 
Specifically the ie prescription of the propagators on the right-hand side comes with an 
opposite sign. 



go 



§ 

o 
o 



n gluon insertions 



o 
o 
o 
o 



m eluon insertions 



Figure 1. Cutting diagram contributing to A mn . The dashed line is the cut. The continuous line 
across the cut in the initial state and in the final state represents a collinear parton. 



(2) The probability P(k±) is related to d 2 A m n/dk 2 L by the equations: 
' d 2 A 

Li< S^-Tri 



P(ki 



-y 

At ^ 

mn 

At^\ dk 2 

mn 



dk 2 
d 2 A 



for k ± ^0 



At ^ 



d 2 k' ± 
(2nf 



d 2 A 

dk' 2 



(3.5) 



for k i 







where A mn stands for A mn averaged over the initial-state polarizations and colors, the 
brackets, (•■•), denote a field average (e.g. a thermal field average) and At is the emission 
time. 6 Eventually we take the limit L — > oo that corresponds to an infinite medium. 

(3) Finally, the jet quenching parameter is obtained by integrating over the probability 
P(k±) as in eq. (3.1). Since in our computation we are considering contributions relevant 



As L provides the length traveled by the parton in the medium, At provides the time. In particular, 
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for a transverse momentum broadening of the parton of order QX or softer, the integral in 
the transverse momentum cannot exceed a cut-off q max 1 which is QX < q nmx <^ Q. 

3.1 Jet broadening in covariant gauge 

The calculation in SCET and in a covariant gauge of P(k±) has been done in [18]. Here 
we sketch a partially different derivation. 

In a covariant gauge, the diagram of fig. 1 involves only vertices of the type (2.6) that 
couple collinear quarks with A + fields of Glauber or soft gluons. Let us first consider the 
case of Glauber gluons only; the diagram is proportional to 

/ II ■ • • 2Qfc+-^ + fc M+te - qi)i 2 Qq t-% ± + ie^ A+{qi ~ q0) ^ M ' 

(3.6) 

where q{ are the collinear quark momenta after each scattering. The Dirac spinor £n(qo) 
satisfies $£n(<Zo) = and is normalized as £^(<7o) £n(?o) = V%Q- As discussed in section 2.3, 
see eq. (2.7), Glauber fields in momentum space do not depend on qf at lowest order in A. 
Hence, the integration over qf yields 

/ dy + d 2 y ± Y[dy^ ■ ■ ■ 6(y^ - yz)A(y + , y% , y±) 9(y^ - y^)A(y + , yf, y±)£,n(qo) , (3.7) 
J i 

where the Glauber fields are now expressed in position space and the functions come from 
the ie prescription in the collinear quark propagators. The same result also holds when 
considering the case of a collinear quark interacting at some point with a soft gluon. For 
illustration, suppose that this happens at the second gluon interaction; the corresponding 
diagram would then contain a term proportional to 

II (2# - ^7-^ + fe - gi¥ 2Qqt l _ ^ + ie U\ qi - ftfrfrfo) , (3.8) 

where A + {q2 — qi) is now a soft field, while A + (q± — qo) is a Glauber field. For qf scales 
as a soft momentum and qf as a collinear one, we may neglect qf ~ QX 2 with respect to 
qf ~ QX in the soft field. Hence, integrating over qf yields again eq. (3.7). 
Equation (3.7) is just a term in the expansion of the Wilson line 



W[y + ,y ± ] = P exp 



L/V2 

ig I dy~A + (y + ,y~,y ± ) 

L/V2 



(3.9) 



with the 6 functions in (3.7) ensuring the path ordering. The distance \^2L is the distance 
that the parton travels along the light cone when moving through the medium in a box 
with a side length L. 



we have 

Af+L 

'~I7T , , At + L 

cte T = =— , 

. **+ L V2 

which is about At/s/2 for At 3> L. The probability P{k±) is then related to an amplitude square normalized 
with respect to the number of particles propagating through the medium, i.e. At/L. 



L 
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From eq. (3.5) it then follows that at leading order in the power counting, and in a 
covariant gauge, the probability P(k±) is given by 

P(k±) = J d 2 x ± e ik ^ ±- (Tr {W*[0,X ± ]W[0,0]}) , (3.10) 

where N c = 3 is the number of colors. The expression holds for the propagation of a 
collinear quark, in which case the gluons in the Wilson lines are in the fundamental rep- 
resentation; the expression for the case of the propagation of a collinear gluon is similar 
and follows from replacing 1/N C by 1/(N 2 — 1) and the gluons in the fundamental repre- 
sentation by gluons in the adjoint one. The trace in (3.10) has to be understood as a trace 
over color matrices. Fields and color matrices are both path-ordered, which, in the plane 
x + = 0, also means that the fields in W[0, 0] are time ordered, while the fields in W^fO, x±] 
are anti-time ordered. 

Although we agree with the result of [18], there are some differences worth emphasizing. 
With respect to the derivation of the result, the step from (3.6) or (3.8) to (3.7) relies on 
the observation that Glauber gluons do not depend, at lowest order, on the momentum 
component along the light-cone direction n. In [18], instead, eq. (3.7) follows from requiring 
QL < 1/A 2 , a request that appears to us neither necessary, nor consistent with integrating 
over energy modes whose momentum components scale like QX 2 . Because in eq. (3.9) we 

I-L/V2 

have chosen a symmetric integration region in y~, taking L — > oo modifies / dy~ 

J-L/V2 

/oo 
dy~ . Note that this limit is not allowed if L is bound by an infrared cut-off. 
-oo 

In the limit L — > oo, the Wilson lines in (3.10) are shown in fig. 2. With respect to the 
interpretation of the result, we note that if soft gluons contribute, their contribution would 
be encoded in the same average of Wilson lines that encodes the contribution from Glauber 
gluons. This observation may be of relevance for perturbative calculations of (3.10) beyond 
leading order. 

(0,-oo, x±) < (0,oo, x±) 



(0,-oo,0) * (0,oo,0) 

Figure 2. The Wilson lines in the right-hand side of eq. (3.10): the upper line corresponds to the 
operator W' [0, x±], the lower line to the operator W[0,0], both taken in the limit L —> oo. 

Finally, we mention that, while W is formally similar to the usual collinear Wilson 
line of SCET, Wn, it has a different origin and the gluon field A + is not collinear. Here, 
W comes from resumming infinite Glauber gluons that interact with a collinear parton 
through a single vertex in the Lagrangian. This leads to a Wilson line integrated along 
the same direction as the momentum of the initial high-energy parton. There, in contrast, 
the gauge invariant building block of SCET [39], n, comes from resumming an infinite 
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amount of n-collinear gluons attached to n-collinear quarks. This leads to a Wilson line 
integrated along the opposite light-cone direction. 

Equation (3.10) was first derived in [17, 25] within different approaches. Clearly the 
equation is not valid in light-cone gauge: if A + (x) = 0, the Wilson line (3.9) becomes equal 
to one and the jet quenching parameter vanishes. This is not surprising since, as discussed 
in section 2.5, new interaction terms show up at lowest order in that gauge and should be 
accounted for. 

3.2 Jet broadening in light-cone gauge 

In this section, we compute the broadening of the transverse momentum of a quark in 
the light-cone gauge A + (x) = 0. This will lead to a gauge-invariant generalization of 
eq. (3.10). We will organize the calculation in an expansion in powers of the background 
(Glauber) gluon fields in the medium, eventually resumming the expansion to all orders. 
An expansion in powers of the gluon fields is sometimes called opacity expansion. First, 
we compute the broadening at leading order in the opacity expansion, i.e. we consider the 
interaction of the collinear quark with one single Glauber gluon. Then we generalize the 
result to an arbitrary order. 

3.2.1 Leading order in the opacity expansion 

We calculate the amplitude An, i.e. the amplitude describing the propagation of a collinear 
quark in a medium of Glauber gluons at leading order in the opacity expansion. The 
corresponding diagram is shown in fig. 3. 



Qo ! k Qo 



> s 


! — § 


— ► — 




: ° 





Figure 3. Leading order diagram in the opacity expansion. 

Computing the diagram of fig. 3 is like computing the imaginary part of the forward 
scattering amplitude of a quark with momentum go = (0,Q,0) on a background field. At 
the cut the quark acquires a transverse momentum k±. Note that, for the calculation is 
performed in light-cone gauge, the quark couples to the transverse component of the gluon 
field and the quark-gluon vertex is at leading order given by eq. (2.15). Gluons can be 
considered as background fields at this stage of the computation. 

Following the definitions given at the beginning of section 3, the contribution of the 
diagram in fig. 3 to d 2 An/d 2 k± reads 

xfn(go) — if- ^ iin{qo)- (3.11) 
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Simplifying the Dirac matrices and writing the gluon fields in position space we find 



d 2 A u 



1 



dk+ dk~ 



2irQ5(2Qk + -kl)e ik < y - y,) 



d 2 k ± * " x± J (2vr)2 

x U%) gAt l {y + \ y-', y±') gAf j (y+, y- } y±)iUio) • (3.12) 



The momentum integrals yield 



dk + dk~ 



(2vr) 2 ---v--- ■w- 



2nQ5(2Qk- ] 



(3.13) 



where the prescription ensures that the function vanishes exponentially in the limit 
|x~| — > oo. Hence, we obtain 



d 2 A 



ii 



d 2 k ± 4V2Q 3 L3 



l — J dy + d 2 y ± d 2 y' ± e- lk ^y^ k\k? L 



dy-'gAT t (y + ,y- ,y±)e- l ^ y 



oo 
oo 



dy gA ± J (y + ,y ,y±)e 2 « 



-Ay -y 'I 



itn{qv) • (3.14) 



Using eqs. (2.9) and (2.12), we may perform the integration over y and y ' obtaining 



d 2 A u 



1 



dy + d 2 y ± d 2 y' ± e - ifc x-(2/W x ) 



d 2 k± V2QL 3 

x Cn(qo) [ g( l> + (y + ,y' ± )-9r(y + ,y±)] [g<P + (y + ,y±)-g(l>-(y + ,y±)] iUqo). (3.15) 

3.2.2 Arbitrary order in the opacity expansion 

The main difficulty in generalizing the previous result to an arbitrary order comes from 
the fact that the interaction between collinear quarks and Glauber gluons can be mediated 
in light-cone gauge by 1-gluon or 2- gluon vertices. The Feynman rule for the leading-order 
1-gluon vertex is given in eq. (2.15), whereas that one for the leading-order 2-gluon vertex 
is in (2.16). The most general diagram contributing to d 2 A nrn /d 2 k± is of the type shown 
in fig. 4. 




Figure 4. Arbitrary order diagram in the opacity expansion in light-cone gauge. Any number of 
1- or 2-gluon vertices can appear on both sides of the cut. 

Let us first consider Feynman diagrams with n gluons on the left of the cut (i.e. the 
part with the ingoing on-shell quark). These are all the Feynman diagrams made of an 
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incoming collinear quark with momentum (0,Q,0), n gluons attached to the quark line, 
and an outgoing collinear quark with momentum (k 2 l /(2Q), k~ , k±), k~ being the sum of 
Q and a residual momentum of order QX or smaller. If n = 1, there is only one possible 
Feynman diagram, the one shown in fig. 3 and computed in the previous section. If n = 2, 
there are two possible diagrams: one made of two 1-gluon vertices and one made of a 2- 
gluon vertex. If n = 3, the possible diagrams are three, if n = 4 the possible diagrams are 
five and so on. In fact, the number of diagrams with n gluons attached to the quark line 
is F n+ i, where F n is the Fibonacci number, implying that the number of diagrams grows 
exponentially for large n. 

An important observation is that, due to the ie prescription of the quark propagator, 
the fields on the left of the cut are time ordered. This means that fields evaluated at 
y~ = — oo (the 4>~ fields) are on the right and fields evaluated at y~ = oo (the 4> + fields) 
are on the left in the expression of the amplitude. We call this amplitude G n (k~ , k±), i.e. 
the sum of all diagrams with n gluons attached to the quark line on the left of the cut. 
Hence, G n {k~ ,k±) can be written as 



G n (k ,k±) 



£ 

3=0 



iQ 



2Qq+ -q{+ie 



(3.16) 



where collects all insertions of (fi + fields and G n all insertions of <j) fields. Once convo- 
luted with the cut and the amplitude on the right of the cut, G n (k±) provides d 2 A nm /d 2 k±. 



First, we calculate G n ; it fulfills the recursion relation 



Gn(q)=l ^U«0x 



+ /(W G -^ X 




(3.17) 



where the Feynman graphs have to be understood as the product of a (momentum space) 
propagator (on the left) and a 1-gluon or a 2-gluon vertex. The solution of the equation is 



G-(q) = / dy+ dy~ 9(-y-yb--®y ++i « + y- f n (y+,q ± )i, 



with 



fn(y + ,q±) 



2Qn\ 



ql / d 2 y ± e- i ^P{[ig^(y + ,y ± )] n ) 



(3.18) 



(3.19) 



The path-ordering P refers to the fields A± that appear in the definition of <fi (see eq. (2.13)): 

/0 rs n -i 
d Sl ...l ds n l ± ■ A ± {x + ,±oo,x ± +l ± 8l) 
-oo J —oo 



x . . . Z_l • A±(x + ,±oo,x± + l±s n ) . (3.20) 



- 15 - 



In particular, it holds that 

V l ± P{[^(x + ,x ± )] n ) =n(Vl0 ± (x + ,x ± ))p([^» ± (x + ,x ± )] n - 1 ) . (3.21) 

The proof of eqs. (3.18) and (3.19) is given in Appendix A. 

The 
relation 

A' 1 



The function G^ can be calculated in a similar manner. It satisfies the recursion 



G+(k-,k ± ,q) 



(2vr) 4 g 



xG+.^AT,^') 



J (2^4 ' o& ' xGt^-fci/), (3.22) 




where the Feynman graphs have to be understood now as the product of a 1-gluon or a 
2-gluon vertex and a (momentum space) propagator (on the right). The solution of the 
equation is 



G+(k-,k±,q) = J dy + dy~e (k 9 )v++ \^ q+ ) V 9(y~) g n (y + ,q ± ,k ± )i, (3.23) 

with 

9n(y + ,q±,k ± ) = ^(kl-ql) f d 2 y ± e-^-^P ([-ig<l ) +(y + ,y ± )] n ) . (3.24) 

The symbol P stands for the anti-path ordering operator and refers to the fields A± that 
appear in the definition of (ft: 

/0 fSn-l 
d Sl ... ds n l ± -A ± (x + ,±oo,x ± + l ± s n ) 

-oo J —oo 

x . . J_l • A±(x + ,±oo,x± + /j_si) . (3.25) 

In this case, it holds that 

V\P{[ ( ft ± (x + ,x ± )] n )=nP^ ± (x + ,x ± )] n - 1 ) {V\(ft ± (x + ,x ± )) . (3.26) 
Substituting the obtained expressions for G~ and into eq. (3.16), we obtain 
G n (k-,k±) = J dy + d 2 y ± S k ~ ~ Q ^ + e ~ ik ^ 

x E -un- i! p ([~ i 9 ( f >+ (y + ^y±)] j ) p ([^r(y + ,y±)] n ' J ) i- (3.27) 

This is the amplitude for n gluons attached to the quark on the left of the cut. The 
amplitude for m gluons attached to the quark on the right of the cut just follows from 
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the Hermitian conjugate of (3.27). Convoluting the two amplitudes and the cut over all 
intermediate momenta, we get d 2 A nm /d 2 k±: 

2vr Q 5(2Qk + - kl) U<lo)Gl(k- , k±) fG n {k~ , fcO&fa,) 



d 2 k ± J (2^) 2 



m ^ 



j'=0 

n 

X 



E -rf - P(M"(y + ,y±)] n " J )^^(go). (3.28) 



From d 2 A nm / 'd 2 k± and eq. (3.5), it follows that the transverse momentum broadening 
probability, P(k±), is given by 

P(k±) = J d 2 x ± e ik ^ x± 

±- (Tr {p ( e -*^-(o.*±)) P (>0 + (o,*x)) p ^ e -i^+(o,o)^ p ( e fc*-(o,o)^ ; (3 29) 



X 



where the trace refers to the color matrices and follows from averaging over the colors of 
the initial state. We have also made use of 

oo n 

P (e"^ + ) P («*♦") = £ £ -^-^ P (l-ig^Y) P ([^~r0 • (3.30) 

n=0 j=0 • / " V ' 

By means of eq. (2.13) we can express P(k±) in terms of the gluon fields A± at light-cone 
infinity: 

P(k±) = J d 2 x ± e lk ^ x ^^(Tr{T\0,-oo,x ± )T(0,oo,x ± )T\0,oo,0)T(0,-oo,0)}^ , 

(3.31) 



where T is the transverse Wilson line [27] 



T(x + , ±oo, x±) = P exp 



->.9 



/•0 

/ ds l± ■ A±(x + , ±oo, x± + /_i_s) 

«/ — oo 



(3.32) 



Note that color matrices and operators are path ordered. The Wilson lines in (3.31) are 
shown in fig. 5. 

3.3 Jet broadening in arbitrary gauge 

In the last section, we have derived in light-cone gauge the probability for a collinear quark 
to gain a certain transverse momentum while travelling through a medium. For the general 
gauge case, one has also to include the interaction between the collinear quark and Glauber 
and soft fields A + , which is encoded in the first term of the SCET Lagrangian (2.2). This 
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(0, — oo, x±) 



(0, oo, x±) 



(0,-oo,0) 



(0, — oo, —oolj 



(0,oo,0) 



(0, oo, -ool±) 



Figure 5. Transverse Wilson lines contributing to the right-hand side of eq. (3.31). For conve- 
nience, in the picture, we have chosen l± \\ x±. 



is done by appropriately extending eq. (3.16) and noticing that in fig. 1 the operators on 
the left of the cut are all time ordered while the operators on the right of the cut are all 
anti-time ordered. As a consequence, on the left of the cut, the transverse contributions 
to always appear to the left and the transverse contributions to G~ always appear to 
the right of the interactions with the A + fields. This allows combining the result obtained 
in section 3.1 in covariant gauge with the result obtained in section 3.2 in light-cone gauge 
to write a fully gauge-invariant expression for P(k±). In the L — > oo limit, it reads 

— oo, x± oo, x± ) T*(0, oo, 0) W[0, 0] T(0, -oo, 0)} \ 

= J rfVe^'^^Trf^^) W(0)}) , (3.33) 

where W has been given in (3.9), T in (3.32), and we have defined W(x±) = Tt(0, 
xW[0,x±] T(0,-oo,x_l) . 

The Wilson lines appearing in eq. (3.33) are shown in fig. 6. Operators and color ma- 
trices are path-ordered along the lines. 7 In the case of non-singular gauges, the transverse 
fields vanish at L — )■ oo, the transverse Wilson lines become one, and the expression of 

7 It has been remarked in [18] that, in the case of equilibrium thermal field averages, a way to handle 
the fact that the fields in the Wilson lines of P(k±) are not time ordered is to express the thermal average 
in the so-called real-time formalism (see e.g. [40-42]). This amounts to modifying the integration path 
along the imaginary-time axis in the partition function to include the real-time axis. More specifically, the 
modification adds to the imaginary-time path a path along the real-time axis at zero imaginary time and a 
parallel path, oppositely oriented, at imaginary time —ie. In our case, the correct ordering of eq. (3.33) would 
be ensured by locating the fields of W(0) on the zero imaginary-time real-axis, and the fields of W (xx) 
on the —ie imaginary-time axis. Fields located at different imaginary times can be treated independently. 
This procedure can be extended to out of equilibrium situations (see e.g. [41, 43]). 
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P{k±) reduces to eq. (3.10). In the case of the light-cone gauge, A + (x) = 0, the Wilson 
lines along the light-cone direction n become one, and the expression of P(k±) reduces to 
eq. (3.31). 

(0, — oo,x ± ) (0, oo,x±) 



(0,-oo,0) 



(0, — oo, — ool±) 



W(0) 
— > 



A (0,00,0) 



(0, oo, — ooZjJ 



Figure 6. Wilson lines appearing in eq. (3.33). Horizontal lines are oriented along the n direction: 
they correspond to the Wilson lines defined in (3.9). Vertical lines extend in the transverse plane 
along the direction of the vector l±: they correspond to the Wilson lines defined in (3.32). For 
convenience, in the picture, we have chosen l± \\ x±. 

It is possible to arrange the Wilson lines in eq. (3.33) in several equivalent ways. First, 
we recall that, according to (2.13) and the following discussion, the transverse Wilson lines 
are independent of the path chosen to connect the initial and the final point. This freedom 
allows us to deform the contour of T(0, oo, x±) such that 



T(0, oo, x±) = [x±, 0±]+ T(0, oo, 0) [—ool±,x± — oo/j_]_ 
where we have defined 



[x±,y±}± = P exp 



o 

g I ds (yj 



x±) ■ A±(0, ±oo,x± + (y± - x_l)s) 



(3.34) 



(3.35) 



Because fields at infinite distance in the transverse plane vanish also in light-cone gauge, 
it holds that [—ool±,x± — oo/j_]+ = 1 and we can write 



P( 



d 2 x ± e ik± - XA 



. (Ti {rt(0, -oo, x ± ) W^[0, x±] [x ± ,0}+ W[0, 0] T(0, -oo, 0)} ) . (3.36) 

Also the transverse Wilson lines at the very left and right of (3.36) combine in a similar 
way, but because the trace only refers to the color matrices and not to the field operators, 
the argument requires some care. It goes as follows. First, we deform the contour of 
Tt(0,-oo,x_L) into ([xj.,0]! = [0,S|_]_) 



rt(o, — oo, xj_) = T^(0, — oo, 0)[0, x±]- , 



(3.37) 
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then we rewrite the trace in (3.36) as 

Tr 1^(0,-00,0) [0,x±]-W^[Q,x±] [x ± , 0]+ W[0,0] T(0, -oo,0)J . (3.38) 
The iV-th term in the expansion of the two transverse Wilson lines in (3.38) is 




xl ± ■ A a f(0, -oo,l ± s n )T a " .. .l ± ■ A a ^(0, -oo,/ ± si)T ai [...] 

/0 rs N rs n+ 2 

ds N / ds N -i ... ds n+ i 
-OO J~ OO J — OO 



xl ± • A a f (0, -oo, l±s N )T aN ...l ± - A a ^ +1 (0, -oo, l±s n+1 )T 




(3.39) 



where the [. . . ] stands for everything that is in between the transverse Wilson lines; the 
limits of integration ensure the proper path ordering. It is now crucial to note that [. . . ] 
contains only gauge fields whose coordinates project in the transverse plane on the straight 
line connecting to x± , while all gauge fields from the transverse Wilson lines are evaluated 
at points that project outside of that line. Moreover, all fields are evaluated at x + = 0. 
These two observations guarantee that the separation between the fields in [. . . ] and the 
fields from the transverse Wilson lines is spacedike. Hence, the gauge fields from the 
transverse Wilson lines commute with [...]. By also using the cyclicity of the trace we can 
then rewrite (3.39) as 

N t-O rs N rs n+2 rO rO r0 

(ig) N y^(-l) N ~ n / ds N / dsN-i... / ds n+1 / ds n / ds n -i . . . / dsx 

n —Q J— OO J —OO J — OO J— OO J S n J S2 

X ; ± • A a " (0, -CO, l±S n ) ...l ± -A ai (0, -CO, l ± Si) 

x l ± ■ A aN (0,-oo, l±3 N )...l± ■ A an+1 (0, -co, l±s n+1 ) 

xTr{[ ] T ajv jia-n+ifan rpaix (340) 

Since all gauge fields in the second and third line are each separated by a spacedike interval, 
they commute as well, and we find 

N_ r0 rs N rs n+ 2 rO rO rO 

(ig) N y^ j (-l) N ~ n / ds N / ds N -i... / ds n+1 / ds n / ds n _i . . . / dsi 

n — Q J —OO J —OO J —OO J —OO J S n J «2 

x l ± -A ai (0, -oo,Jjlsi)...Z_l •^ ajv (0,-co,/_ L S7v) 
x Tr{[. . .}T aN . . .T ai } , (3.41) 

where only the first line depends on the summation index n. The sum adds up to zero. 
This implies that only the zeroth-order term in the expansion of the transverse Wilson lines 
contributes and that T'(0, — oo,0) cancels with T(0, — oo,0) in (3.38). Equation (3.33) can 
thus be rewritten in the equivalent way 

P(k±)= f d 2 x ± e ik ^ x ^(Ti-{[0,x ± ]_W^[0,x ± ][x ± ,0] + W[0,0]}) . (3.42) 
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The advantage of eq. (3.42) is that it is explicitly independent of the choice for the vector 
l± and of the starting point of the integration path in (2.13). The Wilson lines contributing 
to the right-hand side of (3.42) are shown in fig. 7. Note that, while the fields evaluated at 
(0,oo,0) in [rrj_,0] + and W[0, 0] are contiguous, this is not the case for the fields evaluated 
at (0, — oo,0) in [0, x±]- and W[0, 0], for they are separated from each other by fields 
located at light-like distance. In fig. 7, this is signaled by the small gap at (0, — oo,0). 

(0,-oo, xx) < (0,oo, xx) 



(0,-oo,0) 



(0,oo,0) 



Figure 7. Wilson lines contributing to the right-hand side of eq. (3.42). Longitudinal lines are like 
in fig. 6; vertical lines stand for the Wilson lines defined in (3.35). 

We can now prove the gauge invar iance of P(k±). A convenient expression to use is 
eq. (3.36). Under a gauge transformation £l(x), P(k±) transforms as 

— oo, — oo/j_)T^(0, — oo, xj_) W^[0, x_x] 
x[xj.,0] + W[0,0]r(0,-oo,0)n t (0,-c»,-ooZj.)}^ , (3.43) 

where we have put to zero fields evaluated at infinite distance in the transverse plane. The 
cancellation of the gauge transformations to the very left and right follows from noticing 
that they are evaluated at (0, — oo, — oo/j_) and, therefore, commute with all the fields in 
the Wilson lines, since they are separated by space-like intervals. The conclusion is that 
P(k±) defined via eqs. (3.33), (3.36) or (3.42) is gauge invariant. 

Finally, we observe that, although the obtained gauge invariant expression of P(k±) 
reflects expectations on the gauge invariant completion of eq. (3.10) with transverse Wilson 
lines at x~ = ±oo (in the large L limit), the result cannot be found in this form in the 
literature. For instance, an expression similar to (3.33) was found in [25] by extending 
an analogous study on SIDIS to jet quenching. The expression in [25] contains, however, 
only one transverse Wilson line at x~ = — oo but none at x~ = oo. The difference with 
our result might be traced back to the regularization of the singularity in light-cone gauge. 
When choosing an asymmetric prescription, the perpendicular component of the gauge 
field can vanish at either +oo~ or — oo - . In such a case, one of the transverse Wilson lines 
becomes one. 
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3.4 Jet quenching parameter q 

We show now how the obtained expression for P(k±) translates into the jet quenching 
parameter q. According to (3.1) and (3.42), q can be written as 

1 f d 2 k\ 



L 3 J (2vr) 2 



d x±d y±e 



ikx-(xx-y±) 



^V x± V y± (Tr [[y ± ,x ± }.W^[0,x ± }[x ± ,y ± ] + W[0,y ± }}) . (3.44) 



By explicitly evaluating the derivatives acting on the Wilson lines, we find that Glauber 
and soft gluons contribute to the propagation of a collinear quark in a medium along a 
length L — > oo by 

q = V2 [ ^% d 2 x ± dx~e ik ^ 
J (2vr) 2 

x i- ^Tvj [0, x x ]-Ul x [x~ , -oo] gF+ l (0, x~,x ± ) U* x [oo, aT] 

x [x ± ,0]+U 0± [oo, 0] gF+t (0, 0, 0) U 0± [0, -oo] | J , (3.45) 
where F^ 1 = n ■ dA\ — V\A + + ig[A + , A\}. We have also defined the Wilson line 



U x± [x ,y ] = P exp 



i a 



[ dz ^+(0,2! ,x±) 



(3.46) 



which is such that l/l ± [x~ ,y~] = U x± [y~ ,x~] and U x± [L/y/2, — L/y/2] = W[0,x±]. Equa- 
tion (3.45) agrees at leading order in the opacity expansion with a similar expression that 
can be found in [31]. 

The integral over k± in (3.45) has an ultraviolet cut-off, g max , which is of the order of 
QX, the size of the transverse momentum broadening that we are considering. If this cut- 
off could be set to infinity, which may happen in dimensional regularization if the integral 
involves only transverse-momentum regions of order QX or smaller, then the integral in k±_ 
leads to a delta function in the transverse coordinate that squeezes the contour of fig. 7 on 
one line. Under this condition q can be written as 

q = V2 f dx- ir (Tr {[/ 0± [-oo, x~] 5 F+ 4 (0, x~ ,0)U 0± [x~ , 0] 5 F+''(0, 0, 0) C/ 0± [0, -oo]}> . 

(3.47) 

Note, however, that the cut-off cannot be set to infinity in lattice calculations, neither it 
is usually set to infinity in perturbative calculations [10, 32]. 

Finally, when doing perturbative calculations, it is useful to write the medium average 
as an exponential 

■^(Ti^[0,x ± }-W j [0,x ± ][x ± ,0] + W[0,0]^ =e c{x ^ )L , for L -> oo . (3.48) 
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If one only considers the first order in the expansion of the exponential, e c<yX± -^ L ~ 1 + 



C(x±)L + . . . , the Fourier transform of the quantity C(x±) is just the differential rate for 
elastic collisions of a quark with particles in the medium. From (3.1) and (3.42), we thus 
have 



where we have explicitly written the cut-off on the transverse momentum. The left-hand 
side of eq. (3.48) resembles very much a Wilson loop with a transverse extension x± stretch- 
ing along the light-cone coordinate x~ , which, in the x + = plane, is proportional to time, 
whereas — C(x±) resembles the corresponding static energy. The analogy naturally leads 
to the use of lattice data from the static Wilson loop to determine q. However, this re- 
quires some care, for the fields in the usual Wilson loop are time ordered [44], while in 
the left-hand side of eq. (3.48) they are path ordered. We will discuss in the next section, 
in a special case, how this limitation may be circumvented and how we may use already 
existing lattice data to gain information on q. 

4 Application: contribution from the scale g 2 T to q 

As an example for an application of the gauge invariant formulation provided by eqs. (3.48) 
and (3.49), we consider the special case of a jet propagating in a weakly-coupled quark- 
gluon plasma at equilibrium. Because the plasma is weakly-coupled, it is characterized not 
only by the temperature, but also by a hierarchy of other energy scales: the Debye mass, 
gT, and the magnetic mass, g 2 T. While one can compute perturbatively the contributions 
to q coming from the scales T and gT, this is in general not possible for the scale g 2 T, 
not even in the limit g — > [45]. Instead, a way to compute the contributions coming 
from the magnetic mass to the jet quenching parameter is using lattice gauge theories. 
Following [33] , we will argue that the leading contribution coming from the magnetic mass 
is encoded in the static energy of an SU(3) Yang-Mills gauge theory in three dimensions, 
and that it can be extracted from available lattice data. 

We proceed as follows. First, we introduce a cut-off q* to separate contributions 
coming from the momentum region g 2 T from contributions coming from higher-energy 
scales: g max > gT » q* > g 2 T. Then we observe that for momenta below the cut-off and 
up to corrections of relative order g 2 , the ordering of the Wilson lines in the expression 
for q does not matter. The reason can be understood in the real-time formalism [41, 42] 
as due to the fact that, at momenta lower than T, the so-called symmetric propagator is 
the dominant contribution in all two gluon-field correlators, which is a consequence of the 
Bose-enhancement factor. 

At this point we make use of the analysis in [32]. There, it was shown that correlators 
supported on space-like and light-like surfaces, like the left-hand side of eq. (3.48), may 
be analytically continued from Minkowski to Euclidean space-time up to corrections of 
relative order g 2 . After analytical continuation, we furthermore take advantage of the 
hierarchy of thermal scales by systematically integrating them out along the program first 
devised in [46]. Integrating out the temperature in the left-hand side of eq. (3.48) leads 




(3.49) 
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to a thermal field average in a three-dimensional EFT called electrostatic QCD (EQCD). 
The only degrees of freedom of EQCD are the zero modes of the gluon fields, whereas all 
fermionic degrees of freedom and higher modes of the gluon fields have been integrated out. 
Integrating out the Debye mass leads to a thermal average in an EFT called magnetostatic 
QCD (MQCD); MQCD is a three-dimensional theory whose only degrees of freedom are 
the components A , A 2 and A 3 of the gluonic field. In fact, at leading order in the coupling, 
MQCD is exactly an SU(3) Yang-Mills gauge theory in three Euclidean dimensions with 
coupling g 2 D = g 2 T. It is precisely because the coupling provides the only dynamical scale 
of the theory that every Feynman diagram contributes to the same order, and, therefore, 
expectation values in MQCD that do not depend on some external larger scale cannot be 
evaluated in perturbation theory. 

In MQCD, the left-hand side of (3.48) can be read by replacing the field A + with 
A 3 /\/2, since the field A has been integrated out at the energy scale g 2 T. Then, fol- 
lowing [33], we are in the position to relate the contribution to C(x±), coming from the 
magnetic mass, with the static energy in three-dimensional SU(3) gauge theory, V(x±): 

C(x ± ) =-V(x ± ). (4.1) 

9 T 

The static energy V(x±) is a non-perturbative quantity at distances of order l/(g 2 T) that 
has been calculated on the lattice [47]. Note that the identification (4.1) is possible only 
among gauge-invariant quantities, like C(x±) provided by (3.48) is. In terms of V(x±), 
eq. (3.49) can be rewritten as 

^-(^f (4,, 

where Jo (A) is the zeroth-order Bessel function. As suggested in [33], the lattice data 
for (the derivative of) V(x±) can be taken from [47] noticing that a constant shift in the 
potential would not contribute to (4.2). The potential has a short-range tail that can be 
deduced from [48]; a study hat accounts for some features of the short-range potential is 
given in [33]. According to [47], the long-range tail of the potential behaves like 

V(r) = -(a--b^ + ...) , (4.3) 
r \ r r ) 

for r > ro ~ 2.2 /g 2 D . The coefficient a/r 2 with a ~ 1.5 is the string tension and b = 
7r/24 ~ 0.13 is the so-called Ltischer term. By substituting (4.3) into (4.2), we obtain for q 

q* (q*) 3 

the expression a — tt + b This expression, if interpreted as an expansion in o*ro, 

r 2 3 

is convergent for q*r$ < 1, reflecting the fact that it accounts for the long-range part of the 
potential only. A complete study, which is beyond the scope of the present work, would 
require the matching of the long-range part of the potential with a suitable short-range 
part. 

An analysis along the above lines provides just the contribution from the region k± ~ 
g 2 T that enters the perturbative computation of q at NNLO, i.e. at order g e T 3 . In 
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order to have a complete NNLO result, one would also need the contributions from the 
regions k± ~ gT and k± ~ T. What we would like to stress here, however, is that even 
in perturbation theory there are contributions that have to be computed using lattice 
techniques. Such contributions require a fully gauge-invariant definition of q, like the one 
derived in this work. At present, q is known up to NLO in the g expansion [32] (see also 
[49] for a leading-order analysis of eq. (3.10) in different transverse momentum regimes), 
so non-perturbative physics will indeed be needed starting from the next order. 

5 Conclusions 

We have derived a gauge invariant definition of the jet quenching parameter q under the 
assumptions that the medium is very large (L — > oo) and that the jet energy, Q, is much 
larger than any other energy scale of the medium (e.g. Q S> T). The existence of very 
different energy scales allows for the construction of an EFT, namely SCET supplemented 
with Glauber gluons, that describes the transverse momentum broadening of a highly 
energetic particle due to the medium. The effective theory is organized in a systematic 
expansion in A, which is the small parameter associated with the ratio of the low energy 
scales of the medium and Q. The specific power counting of the SCET Lagrangian (2.2) 
depends on the gauge: some terms are enhanced in light-cone gauge with respect to a 
covariant gauge. Hence, additional vertices need to be considered in the general gauge case 
with respect to the simpler covariant gauge case. By a direct calculation of all diagrams 
containing these additional vertices at lowest order in A, we have found that q is related 
to the medium average of some previously known longitudinal Wilson lines as well as 
some transverse Wilson lines. The main result for the transverse momentum broadening 
probability is given in eq. (3.33). The corresponding Wilson lines are shown in fig. 6. 

With respect to other expressions that can be found in the literature, eq. (3.33) appears 
to hold for all regularizations of the light-cone singularity in light-cone gauge. Moreover, 
it follows from an explicit resummation of Feynman diagrams and not from a heuristic 
extension of the covariant gauge result. The fields in the Wilson lines of (3.33) are path 
ordered, which means that fields supported in one transverse plane are time ordered while 
fields supported in the other transverse plane are anti-time ordered. This leads to some 
subtleties when explicitly proving the gauge invariance of the expression; a proof of gauge 
invariance can be found at the end of section 3.3. Equivalent formulations of eq. (3.33) are 
given by eq. (3.36) and eq. (3.42). The Wilson lines contributing to the latter are shown 
in fig. 7. 

The fully gauge-invariant expression not only allows for computations in any gauge, 
but also opens the way for the use of lattice data in the evaluation of q, as suggested in 
[30, 32, 33]. A particularly suitable expression is provided by the transverse momentum 
broadening probability given in eq. (3.42), which translates into the expression of q given 
in eqs. (3.48) and (3.49). As an example of using lattice data, we have discussed the 
contribution from the momentum region k± ~ g 2 T to q in a weakly-coupled quark-gluon 
plasma. 
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The calculation of the jet broadening presented here includes the effect of Glauber 
and soft gluons. There are, however, other modes whose contributions may be relevant 
and that will need to be considered before comparing with data. The addition of these 
new modes, in particular collinear modes [35], may be systematically accounted for in the 
SCET framework. Also the calculation of q in a weakly-coupled quark-gluon plasma may 
be improved by including the remaining (perturbatively calculable) NNLO terms [50]. 
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A Computation of G^(q) 

In this appendix, we derive eqs. (3.18) and (3.19). The proof of eqs. (3.23) and (3.24) is 
analogous and, for this reason, will not be detailed here. 

We will proceed as follows. First, we prove eqs. (3.18) and (3.19) for the cases n = 1 
and n = 2. After that, we will show that if the relations (3.18) and (3.19) are fulfilled 
for n — 2 and n — 1, then they are fulfilled for n, which proves eqs. (3.18) and (3.19) by 
induction. 

A.l Computation of G^(q) 

The computation of Gi(q) is straightforward; by just writing down the Feynman rule we 
obtain 

n-( \ iq± " 9 A i n (Q - go) / / a -i \ 

G 1 (q) = r ■ ( Al ) 

Expressing the gluon field in position space and using eq. (2.12), we can then write 

GT(q) = ~^qJ d y + dy~e(-y~)e^--®y ++t o + y- q\ #"(y + ,^) ? /. (A.2) 



Comparing with eq. (3.18), eq. (A.2) implies 



2Q 



g4> (y + ,q±), 



(A.3) 



which agrees with (3.19) for the case n = 1. 

A.2 Computation of G^iq) 

From eq. (3.17), it follows that G^iq) satisfies 



go 




+ 



(2vr) 4 



GT(gi)x 



(A.4) 
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Hence, G 2 (g) gets a 2-gluon vertex contribution, namely 




(A.5) 



and a 1-gluon vertex contribution, 



G 2b (q) = J 



Al 
(2vr) 4 



G 1 (q x ) x 



(A.6) 



Computing the two, we obtain 



G- 2a {q) 



I d 2 q 1± 



2QJ (2vr) 2 



/ d 2 q l± 



2QJ (2tt)< 



dy + dy-0{-y-) e ^ q ~ Q ^ ++iq+y 

x {q± - qi±) ■ qi± g^~(y + ,q± - qi±) g4>~(y + \qi±)i, 

(A.7) 

dy + dy-0(-y-) e *(9~-0)i/ + +<9 + »~ 

x (?i-9i±) 94>~{y + ,q± - qi±) g<l>~(y + ,qi±)$, 

(A.8) 



and finally G 2 (g), which is of the form (3.18), but with f 2 given by 



f2(y + ,q±) = ^ 



I d 2 qi± 



2QJ (2vr) 2 



[(q±- qi±) ■ qi± - (q±-qf±)] 

x g<p~{y + ,q± - q 1± ) g4>~ (y + , gi ±) . (A.9) 



In order to show that eq. (A.9) is equivalent to (3.19), we express (ft in position space and 
use eq. (3.21): 

f2(y + ,q±) = Trk / l^lo gi • (gi± - q±) g^~{y + ,q± - gn) #"(y + ,gi±) 



2Q7 (2vr) 2 



2Q 



d 2 y± d 2 x± 



d 2 gi± 
(2k) 2 



q± ■ -iV±e 



x g4> (y + ,y±) g<P (y + ,%±) 



2Q 



d 2 y±e - iC >^q ± • {-iV±gr(y + ,y±)) #~(y + ,y±) 



— q\ I 4e^P 



(y + ,y±)] 2 ) , 



(A.10) 



which agrees with (3.19) for the case n = 2. 
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A. 3 Proof by induction 

In order to complete the proof by induction, we need to show that if G~_ 1 and G~_ 2 fulfill 
eqs. (3.18) and (3.19), then G~, defined through eq. (3.17), also fulfills them. Specifically, 
from eq. (3.17) it follows that G~ is of the form (3.18) with f n given by 

d 2 q n -i± gj - gn-i. 
2vr) 2 ql 

d 2 q n - 1± d 2 q n - 2 ± (g± - q n -i±) ■ (gn-i± - q n -2±) 



fn(y + ,q±)=i —F7TT2 2 9<f> {v , Q± ~ 9n-l ±) fn-i (y , gn-l ± ) 



(2vr)2 (2vr)2 q \_ 2 

x g<P~(y + ,q±-qn-i±) g<P~(y + ,q n -i±-q n -2±) fn-2(y + ,qn-2±)- (A.n) 

Using the expressions of / n -i and / n -2 given in eq. (3.19), we get 

fn(y + ,q±) = - y d d 2 y±n (gi - gLi±) 9^~{y + ,q± - q n -i±) 



* /" d 2 gn-l± ^ 2 gn-2± ,2 / i\ / \/ ^ 

+ 2^! y "w - "w - y± n(n " } (g± " qn ~ i±] {qn ~ i± " gn_2±) 

x g<P~{y + ,q± - q n -i±) g<P~(y + ,q n -i± - q n -2±) 

^p([^-(y+,y ± )] n - 2 ) . (A.12) 



Finally, by using similar manipulations as in (A. 10), this can be brought into the form of 
eq. (3.19). 
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